Geometric resonances in the magnetoresistance of hexagonal lateral superlattices 
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We have measured magnetoresistance of hexagonal lateral superlattices. We observe three types 
of oscillations engendered by periodic potential modulation having hexagonal-lattice symmetry: 
amplitude modulation of the Shubnikov-de Haas oscillations, commensurability oscillations, and 
the geometric resonances of open orbits generated by Bragg reflections. The latter two reveal 
the presence of two characteristic periodicities, v3a/2 and a/2, inherent in a hexagonal lattice 
with the lattice constant a. The formation of the hexagonal-superlattice minibands manifested 
by the observation of open orbits marks the first step toward realizing massless Dirac fermions in 
semiconductor 2DEGs. 

PACS numbers: 73.43.Qt, 73.23.-b, 73.21. Cd 



I. INTRODUCTION 

A hexagonal lateral superlattice (HLSL) — a two- 
dimensional electron gas (2DEG) subjected to periodic 
potential modulation with hexagonal-lattice symmetry 
— is of interest in two different contexts. First, it is en- 
visaged as a route to artificially generate massless Dirac 
fermions (MDF) at the corners of the superlattice Bril- 
louin zone.— ~— Second, it is expected to stabilize^ the frag- 
ile "bubble phase" (the hexagonal-lattice arrangement of 
two- or three-electron clusters) in the quantum Hall sys- 
tem, which has been theoretically predicted to be the 
ground state at the 1/4 or 3/4 fillings of the third or 
higher Landau levels^— Analogous stabilization is re- 
ported for the stripe phase at the half fillings, using 
one-dimensional (ID) lateral superlattices . 10 ! 11 As an ini- 
tial step toward pursuing these intriguing possibilities, 
we study, in the present work, low-field magnetoresis- 
tance of lateral superlattices with a weak (a few per- 
cent of the Fermi energy E-p) hexagonal-lattice potential 
modulation, fabricated from conventional GaAs/AlGaAs 
2DEGs. 

Numerous studies have been devoted to lateral super- 
lattices having lDi2r— and two-dimensional (2D) square- 
±±JH or rectangular-lattic o 21 ' 22 potential modulations. 
By comparison, magnetotransport of HLSLs remains rel- 
atively unexplored. An early experiment by Fang and 
Stiles^ 3 - exhibited, for a weak modulation amplitude, 
commensurability oscillations (CO) similar to those ob- 
served in ID lateral superlattices^ 2 - but corresponding to 
the periodicity half the lattice constant of the hexag- 
onal lattice. Hexagonal lattices are not uncommon in 
antidots,— ~— which represent the strong limit of the 
modulation amplitudes. As demonstrated in Ref . l23l (and 
also in Ref. for square lattices), however, antidots and 
weakly modulated lateral superlattices display qualita- 
tively distinct behavior. 

In the present paper, we report three variants of os- 
cillatory phenomena in HLSLs qualitatively similar to 
those known in ID lateral superlattices: CO (some- 
times also called as "geometric resonances") due to the 



commensurability between the cyclotron radius and the 
periodicity in the superlattice,— amplitude modulation 
(AM) of the Shubnikov-de Haas oscillations (SdHO) ,22^ 
and an alternative type of geometric resonances, geomet- 
ric resonances of open orbits (GROO) i^2r— Open orbits 
are composed of segments of cyclotron orbits repeatedly 
diffracted by Bragg reflections from the superlattice po- 
tential. Resonances take place when the width of the 
open orbits coincides with the periodicity. In both CO 
and GROO, oscillations are observed as superposition of 
two components originating from the periodicities y3a/2 
and a/2, respectively, with a representing the lattice con- 
stant of the hexagonal lattice. The CO for the latter 
periodicity corresponds to that observed in Ref. [23| men- 
tioned above. We will discuss the amplitude of CO in 
connection with the amplitude of the potential modula- 
tion inferred from the AM of SdHO. 

One of the major motivations in the exploration of su- 
perlattices is to artificially design and generate a band 
structure (miniband) that possesses the length and en- 
ergy scale different from that in natural crystals^ with 
the generation of MDF in HLSL being one example. In 
lateral superlattices, however, clear evidence of the for- 
mation of the miniband structure has not been observed 
until recently) 38 ' 39 probably due to the technical difficulty 
in fabricating, with minimal disorder, superlattices with 
the period small enough (close to the Fermi wavelength) 
to embrace minibands. Note that both CO and AM of 
SdHO can be traced back to the oscillation with the mag- 
netic field of the width of the Landau bands (Landau lev- 
els broadened by the modulation potential), and do not 
require minibands for their occurrence. By contrast, the 
observation of GROO evinces the formation of miniband 
structure, attesting to the high quality of our superlattice 
samples eligible to seek for artificial Dirac fermions. 

The paper is organized as follows. After describing 
experimental details in Sec. UH experimentally obtained 
magnetoresistance traces exhibiting CO, AM of SdHO, 
and GROO are presented in Sees. UlLAl iHIBl andHTTCl 
respectively. Prospects and necessary improvements to 
be made to realize MDF and to detect the bubble phase 
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FIG. 1. (Color online) Hexagonal lateral superlattice (HLSL) 
sample used in the present study, (a) Hall bar containing 
both the section with hexagonal-lattice potential modulation 
(HLSL) and the section without modulation (plain 2DEG) for 
reference, (b) Schematic drawing of the HLSL section, (c) 
Scanning electron micrograph of the hexagonal lattice of EB- 
resist (black dots) placed on the surface of the 2DEG wafer. 



are discussed in Sees. IIV Al and IIVB1 respectively, fol- 
lowed by concluding remarks in Sec. [V] 



II. EXPERIMENTAL DETAILS 

Schematics of the HLSL samples used in the present 
study are depicted in Fig.[TJ The samples were fabricated 
from a conventional GaAs/AlGaAs 2DEG wafer with the 
heterointerface residing at the depth d = 60 nm from the 
surface. As shown in Fig.QJa), the 2DEG wafer was pat- 
terned into Hall bar with the width 40 /im and contain- 
ing two sets of voltage probes (with the inter-probe dis- 
tance 320 /im) to measure the section with (HLSL) and 
without (plain 2DEG, for reference) the hexagonal-lattice 
modulation. The potential modulation was introduced 
by placing a hexagonal lattice of high-resolution nega- 
tive electron-beam (EB) resist (calixarene) 4 ^ on the sur- 
face (Fig. [Tfb)(c)) and making use of the strain- induced 
piezoelectric effect^ as has been done to prepare ID 
lateral superlatticeSi 10 ! 11 ' 33-36 ' 42 ^ 3 Compared with more 
general methods to introduce potential modulations, e.g., 
by placing metallic gate grids or by shallow etching, the 
simpleness of our approach (only one-step process, the 
EB drawing, is needed to introduce the modulation), 
along with the high spatial resolution of the EB resist 
we employ, allows us to prepare highly ordered lateral 
superlattice samples with minimal damage to the 2DEG. 
In fact, the mobility p = 88 m 2 V _1 s _1 and the electron 
density n e = 3.9 x 10 15 m~ 2 of the 2DEG wafer remained 
virtually unchanged after the fabrication of the HLSL de- 
vices. We prepared HLSL samples with the lattice con- 



stant a — 200 nm and 100 nm. Since the modulation 
amplitude for a = 100 nm was found to be extremely 
small, we mainly present the data taken from a = 200 
nm HLSL in the followings. Note that the modulation 
strength rapidly decreases with decreasing a, roughly as 
cxp(— a/d). 43 Resistivity measurements were carried out 
employing standard low-frequency ac lock-in technique at 
4.2 K for CO and GROO, and at 15 mK, using a dilution 
refrigerator, for SdHO. 



III. EXPERIMENTAL RESULTS 

A. Commensurability oscillations 

Magnetoresistance Ap(B)/po of a HLSL with a = 200 
nm is shown in Fig. [2]Ja). Here, po is the resistivity p(B) 
at B = and Ap(B) = p(B) — po. Oscillatory behav- 
ior is apparent above ^0.1 T (see also the bottom trace 
in Fig. Etc)), which is absent in the trace for the plain 
2DEG. Small oscillations seen above ^0.35 T for both 
traces are the SdHO. To analyze the oscillations, we take 
the second derivative (d 2 /dB 2 )(Ap(B)/ p Q ) numerically 
to extract the oscillatory part, plot it against 1/B (inset 
to Fig. [UJb)), and then perform Fourier transform. The 
Fourier spectrum thus obtained is presented in the main 
panel of Fig. Elb). Three eminent peaks, Ja, /b; an d 
/sdH, are seen: /sdH = n e h/(2e) represents the SdHO, 
while /a and /b coincide with the frequency 2fi,fcp/(ep), 
with = ^/2Tm e the Fermi wave number, of the CO 
corresponding to the periodicity p = \f2>a/2 = pa and 
a/2 = p B , respectively. As depicted in Fig.[3{a), bothpA 
and pb are representative spacings between the lattice 
points contained in the hexagonal lattice. More gener- 
ally, lattice spacings in the hexagonal lattice are given 
by P(h,k) = 2n/\g( h ,k)\, with p A , pb, and p c correspond- 
ing to (h, k) = (1,0), (1, 1), and (2, 1), respectively. Here, 
Q(k,k) — ha\ + ka\ represents a reciprocal lattice vector 
with a\ = (27r/a)(l,-l/v / 3) and a* 2 = (2ir/a)(l,l/V3) 
the primitive reciprocal lattice vectors. A potential mod- 
ulation having the hexagonal-lattice symmetry can gen- 
erally be written as, 



V ( r ) = H y P( h , k) cos{g (htk) ■ r), 



(1) 



h.k 



with r = (x,y). We defined the direction of the current 
in the measurement of the resistance as the x axis (see 
Fig. Hlja)). The summation is taken over sets of inte- 
gers (h, k) that yield independent periodicities in the 2D 
plane. 44 Hexagonal modulation potential often given in 
the literatures (e.g., Eq. (3) in Ref. l23l and Eq. (1) in Ref. 
0, see also Eq. (J7|) below) corresponds to Eq. (JJ) retain- 
ing only the components with the largest spacings pa, 
namely, V p with (h,k) = (1, 0), (1, -1), (0, -1). The 
presence of the two peaks /a and /b in the Fourier spec- 
trum reveals that the observed CO is the superposition 
of two components corresponding to the periodicities pa 
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FIG. 2. (Color online) (a) Magnetoresistance traces at 4.2 
K for the HLSL with a = 200 nm (solid line) and the ad- 
jacent plain 2DEG (dashed line). Dotted curve represents 
slowly-varying background for the HLSL section, obtained by 
polynomial fitting, (b) Inset: Oscillatory part for the HLSL 
section obtained by taking the second derivative with respect 
to B, plotted against jB -1 . Main: Fourier spectrum of the 
oscillatory part shown in the inset, (c) Oscillatory part of 
the magnetoresistance for the HLSL section obtained by sub- 
tracting the slowly-varying background shown in (a) (bot- 
tom), components having the frequency /a and /b extracted 
by the Fourier bandpass filter method (see text for details), 
and the addition of the /a and /b components (top). The 
latter traces are offset by 0.002 for clarity. Vertical ticks for 
the /a and /b components indicate the positions of the n-th 
flat-band conditions given by Eq. ([2}. Dotted curves are the 
fit to Eq. ©. 
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FIG. 3. (Color online) (a) Representative lattice spacings in 
the hexagonal lattice. The periodicities pa and pb account for 
the components /a and /b in Fig. O respectively. Periodici- 
ties equivalent to those shown are found by rotating them by 
±120°. (b) Drift velocity Vdi along the direction 8 from the 
y-axis and «d2, «d3 along the other two equivalent directions. 
9 = 30°, 0°, and arccot(3\/3) — 10.9° for the periodicities pa, 
Pb, and pc, respectively. 



and pb , and that the components corresponding the pe- 
riodicity p B , V Pih h) with (h,k) = (1,1), (2,-1), (1,-2), 
are also present in the modulation potential Eq. ([T]). 

To gain more quantitative information from the CO, 
we separate out individual components following the pre- 
scription described in Ref. [jot first we perform Fourier 
band-pass filter on (d 2 / dB 2 ){Ap{B) / p ) vs. B^ 1 plotted 
in the inset of Fig. [2j employing the window that encom- 
passes the peak /a or /b (shaded areas in Fig. [UJb)); 
(d 2 /dB 2 )(Ap(B)/po) corresponding to a single compo- 
nent thus extracted are replotted against B and then 
numerically integrated by B twice. The traces of sin- 
gle component (Ap(B)/p )co f° r /a and /b restored by 
this procedure are plotted in Fig. [2c). As demonstrated 
in the figure, addition of the /a and /b components re- 
produces well the oscillatory part of Ap(B)/p directly 
obtained by subtracting the slowly varying background 
found by polynomial fitting to the data (plotted with the 
dotted line in Fig. EJa)), apart from the SdHO above 
^0.35 T. In Fig- EKc) , positions of the flat-band condi- 
tions, where the drift velocity vanishes, 
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are indicated by vertical ticks. Here, R c = hkp/(eB) is 
the cyclotron radius. It can be seen for both components 
that the minima take place at the flat-band conditions, 
as is the case with ID lateral superlattices . 12 ' 46 To be 
more precise, the diffusion (band) contribution resulting 
from the drift velocity takes minima, while the colli- 
sional (hopping) contribution due to the modulation of 
the density of states (DOS) takes maxima, at Eq. ([2]) 
in ID lateral superlattices^ with the former being by 
far the dominant contribution in the measurement of the 
resistivity along the principal axis of the modulation. 

We further make an attempt to fit the extracted CO 
curves to a formula representing the diffusion contribu- 
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tion for ID lateral superlattices V(x) = V p cos(2ttx/p), 
in which the damping of the oscillations by scatterings is 
taken into account : 42 ! 48 

(po)co lA ( T p) A (mw-b) p V P Bsm ( 27r p )' 

(3) 

where 

7 2(2^)3/2 { e J { 2 m* J n 3/2' 1 ' 

and T p = pfcFSwc/(47r 2 fcB), with m* the effective mass, 
fee the Boltzmann constant, w c = eB/m* the cyclotron 
angular frequency, and A(X) = X/ sinhA. We use the 
modulation amplitude V p and the effective mobility /i w 
as fitting parameters. As shown by dotted curves in Fig. 
[^c), fairly good fitting is achieved, albeit within rather 
limited magnetic field range: B < 0.20 T for Ja and 
B < 0.37 T for The values of V p and p w obtained by 
the fittings are V p = 0.025 and 0.016 meV, and p, w = 9.1 
and 8.2 m 2 V _1 s _1 for /a and /b (p — Pa — 173 nm and 
Pb = 100 nm), respectively. It has been shown for ID 
lateral superlattices 4 ^ that /z w is close to the single par- 
ticle (or quantum) mobility p s that describes the damp- 
ing oc exp(— ir/p s B) of the SdHO. 50 This is found to be 
roughly the case also for our HLSL; we obtain /i s = 6.2 
m V s from the similar Fourier analysis of SdHO for 
the data shown in Fig.^a) (and also for the data taken at 
T ~ 15 mK, see Fig. [5] dependence of p s on the temper- 
ature was not observed in this temperature range). The 
values of Vp, on the other hand, appear to be too small. 
Much larger values were found for ID lateral superlat- 
tices having similar modulation periods and fabricated 
from the same 2DEG wafer and therefore expected to 
have similar modulation amplitudes: 0.18 and 0.06 meV 
for the periods a = 200 and 90 nm, respectively] 51 ' 52 
The V p obtained here for the HLSL cannot be literally 
taken to represent the modulation amplitude in Eq. ([T]) 
correctly for the following reasons. 

First, it is necessary to note that the observed CO is 
an addition of contributions from three equivalent mod- 
ulations with the same spacing present in the hexagonal 
lattice, rotated by 120° from each other. The drift ve- 
locity Vd(<x V p ) responsible for the diffusion contribution 
is pointed perpendicular to the direction of the modula- 
tion axis. As depicted in Fig. [31(b), the drift velocity is 
directed toward the angle 9 (v dl ) and 9± 120° (v d 2) v d3 ) 
from the y axis, with 9 = 30° and 0° for the modulations 
A and B, respectively. Since p xx ~ Vyy/Vyx (with cr Q( 3 
representing components of the conductivity tensor), CO 
is proportional to the y component of the drift velocity 
squared, Ap C o oc |v d i| 2 cos 2 9 + |v d2 | 2 cos 2 (6> + 120°) + 
|v d3 | 2 cos 2 (6> - 120°). This equals (3/2)u| regardless of 
the angle 9 if we assume |v d i| = |v d2 | = |v d3 | = var^- 
leading to the correction of V p to the factor of y/2/3 
smaller values. 

More importantly, it has been shown that the ampli- 
tudes of CO for 2D lateral superlattices are usually much 



smaller than those for ID lateral superlattices having a 
similar modulation amplitude^ This was initially at- 
tributed to the splitting of the Landau levels into sub- 
levels (Hofstadter spectru m 54 ' 55 ), which suppresses the 
diffusion contribution! 16 ' 18 Later, an alternative explana- 
tion was presented by Grant et a/.— based on the calcula- 
tion of semiclassical trajectories of the electrons showing 
that the drifting motion introduced by the modulation 
in the x direction is suppressed by the modulation in the 
y direction. They showed that the diffusion contribu- 
tion survives, with the amplitude reduced compared to 
ID modulation, only when the modulation is asymmetric 
between x and y directions. For symmetric modulation, 
the diffusion contribution vanishes, leaving only the colli- 
sional contribution having the oscillation phase opposite 
to that of the diffusion contribution. The effect of the 
symmetry between x and y directions was experimentally 
confirmed.— Although these are for a 2D square lattice, 
qualitatively similar mechanism is expected to be opera- 
tive also in the hexagonal lattice, in which modulations 
with differing orientations coexist. Therefore the value of 
V p obtained by fitting Eq. Q to the CO trace is expected 
to be smaller than the modulation amplitude also in the 
hexagonal lattice. This is confirmed by comparing V p 
obtained here to the modulation amplitude inferred from 
the AM of the SdHO, as will be shown in the subsequent 
subsection IIIIBI 



B. Amplitude modulation of Shubnikov-de Haas 
oscillations 

In Fig. 0Ja), we plot Ap/po taken at T = 15 mK. 
It can readily be seen that the SdHO for the HLSL ex- 
hibits modulation in the oscillation amplitude, with the 
amplitude maxima at the flat-band conditions Eq. ([5]) 
for the periodicity pa = v3a/2. The AM is also evi- 
dent in the Fourier spectrum shown in Fig. EJb) taken of 
(d 2 /dB 2 )(Ap/ po) vs. B~ x , which exhibits, in addition 
to the peak representing the principal frequency /s d H of 
the SdHO and its harmonics, side peaks marked with 
/a+ and /a-, with the distance |/a± — /sdH| being equal 
to the frequency of the CO, 2hk-p / '(epA), corresponding 
to the periodicity pa; the distance coincides with /a in 
Fig. (Ub) after the correction of the small difference in 
the electron density n e between different cooling downs. 
The AM attributable to the periodicity p-Q was not clearly 
observed, probably owing to the weakness of the modu- 
lation. Note that, as mentioned earlier, the amplitude 
of the potential modulation rapidly decreases with de- 
creasing periodicity. The absence of the ps component 
indicates that the AM of SdHO is more heavily weighted 
to larger amplitude components of the potential modu- 
lation compared to the CO. 

The AM of SdHO is known to also originate from the 
two mechanisms, the diffusion contribution and the col- 
lisional contribution, with the amplitude minima (max- 
ima) taking place at the flat-band conditions for the for- 
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FIG. 4. (Color online) Magnetoresistance traces at 15 mK for 
the HLSL with a = 200 nm (solid line) and the adjacent plain 
2DEG (dashed line, offset by 1.0 for clarity), showing rapid 
SdHO. Slowly- varying background obtained by averaging the 
upper and lower envelop curves of the Ap/po for the plain 
2DEG is shown by dot-dashed line. (The same background is 
shown for both traces). Vertical ticks indicate the flat-band 
conditions, Eq. ([2]), for periodicities pa and pb- (b) Fourier 
spectra taken of (d 2 /dB 2 )(Ap/p ) vs. B" 1 for the HLSL 
(solid line) and the plain 2DEG (dashed line) sections. 



mer (latter) mechanism.—^ For ID lateral superlattices, 
it has been shown that the collisional contribution dom- 
inates at low magnetic fields (< 0.25 T)^ This appears 
to be also the case in our HLSL, as can be seen in Fig. 
IDJa) exhibiting amplitude maxima at the flat band condi- 
tions. The diffusion contribution in the SdHO in HLSL, 
if any, is expected to be much smaller than in ID lat- 
eral superlattices, by analogy with the case for the CO 
discussed above. 

The collisional contribution to SdHO in ID lateral su- 
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FIG. 5. (Color online) Oscillatory parts of the magnetore- 
sistance shown in Fig. [3] obtained by subtracting the slowly- 
varying background (bottom) and those calculated by Eqs. 
© and © with n e = 3.9 x 10 15 mT 2 , fi B = 6.2 m 2 V -1 s-\ 
T = 15 mK, C — 2.2, pa = 173 nm, and the indicated values 
of V PA for the plain 2DEG (a) and the HLSL (b) sections. 
Calculated traces are offset by 1.0 for clarity. 



perlattices V(x) = V p cos(2irx /p) is described well by. 
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representing the Landau bandwidth, T c = hiu c /(2ir 2 kB), 
C a constant ^2,— and Jo(x) the Bessel function of order 
zero. According to Wang et al.^ bandwidth of the N- 
th Landau level for a 2DEG subjected to the hexagonal 
potential modulation^ 

V(r) = V PA [cos(al ■ r) + cos(a2 • r) + cos^ • r)] (7) 

with 0,3 = — 0.3 is given by jV PA e~ u / 2 LKr(u), where 
u = (2tt/pa) 2 1 2 /2 with I = ^JK/JeB) the magnetic 
length, and Ln(u) represents the Laguerre polynomial. 
The Landau bandwidth at the Fermi energy at a low 
magnetic field (N 3> 1) is thus approximated well by 
simply 9/4 times Eq. (J6)): 



9 1 [pa ( 2tt tt 
Vb = t^pa-W -5- cos — R c - - 
4 R c \pa 4 



(8) 



We therefore make an attempt to analyze AM of SdHO 
shown in Fig. |4ja) using Eqs. §5$ and dSJ. In Fig. [5j we 
compare traces calculated by Eqs. ([5]), ([5} with SdHO ex- 
tracted from the experimental Ap/po shown in Fig.^a) 
by subtracting the slowly-varying background. As can be 
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seen in Fig.[5la), the trace calculated with V PA = repro- 
duces the experimental trace for the plain 2DEG section 
quite well. Figure EJb) shows that AM is barely visible 
if we use the value V PA — 0.021 meV obtained by fitting 
the CO trace to Eq. ([3]) and applying the correction for 
the factor y/2/3 to account for three different directions 
of the drift velocity outlined in Sec. MI Al To reproduce 
experimental AM, a much larger value V PA =0.12 meV 
is required. Note that Eq. ([5} already includes the con- 
tribution from the three equivalent orientations (see Eq. 
(O). Since the collisional contribution is the effect of 
modulated DOS that alters the scattering rate of elec- 
trons, it essentially does not depend on the direction of 
the modulation^ and therefore interference between dif- 
ferent directions as in the diffusion contribution is con- 
sidered to be absent. We thus expect the value of V p 
derived from the analysis of the SdHO presented here to 
represent better the amplitude of the modulation. The 
value is also consistent with the modulation amplitude of 
ID lateral superlattices with similar modulation periods 
mentioned earlier. 



C. Geometric resonances of open orbits 

We find yet another type of oscillations in the 
magnetic-field range below the regime where CO is ob- 
served. The small amplitude oscillations become clearly 
visible after the subtraction of slowly-varying back- 
ground, as illustrated in Fig. [HI Figure [^b) reveals that 
the oscillations observed at B > are reproduced, in- 
cluding minute structures, in the trace taken at B < 0, 
ruling out the possibility that they simply result from 
external noises. Since the measurement was performed 
on a large (40 x 320 /im 2 ) Hall bar at a relatively high 
temperature, T = 4.2 K, the oscillations are unlikely to 
be related to the well-known universal conductance fluc- 
tuations (UCF) i 60 i 61 The patterns of the oscillations do 
not change between separate cooling downs (see Fig. [9] 
below), which is also at variance with UCF. Similar 
small amplitude oscillations have been reported in ID 
lateral superlattices, and interpreted as the geometric 
resonances of open orbits (GROO)i^— The oscillations 
observed here in HLSL can also be interpreted with the 
same mechanism. 

Open orbits are generated by Bragg reflections due to 
the superlattice potential. Figure [T^a) shows Fermi cir- 
cles in the reciprocal space in the extended zone scheme, 
with the open orbits depicted as chains of Fermi circle 
segments (thick curves). As illustrated in Figs.[7fa) and 
(b), we denote the open orbit composed of arcs from j-th 
and fc-th nearest-neighbor Fermi circles as X(j, k), with 
X = A, B for the the open orbits engendered by the pe- 
riodicity pa and pb, respectively. The open orbit X(j, k) 
is generated by the Bragg reflections from j-th and fc-th 
harmonics of the modulation potential with the period- 
icity px- Open orbits in the real space are obtained by 
rotating those in the reciprocal space by 90° and mul- 
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FIG. 6. (Color online) (a) Low-field magnetoresistance traces 
of the HLSL section (a = 200 nm) at 4.2 K for B > and B < 
(offset by l.OxlO -3 ), plotted against the absolute value \B\. 
Slowly-varying backgrounds acquired by polynomial fitting 
are shown by dot-dashed lines, (b) Oscillatory parts for B > 
and B < (offset by 0.2xl0 -3 ), obtained by subtracting 
the slowly- varying background. Downward triangles mark the 
positions for the transverse resonance of the orbit X(j, fc) (X 
= A or B), Eq. ([9]). Vertical ticks indicate the positions for 
the n-th longitudinal resonance, Eq. (|11[1 . 



tiplying by the factor H/(eB), as portrayed in Fig. Etc). 
The width &xo,fc) °f the orbit X(j, k), therefore, decreases 
inversely proportional to B. The transverse resonance 
takes place, leading to the maxima in the conductivity 
(hence in the resistivity as well), at magnetic fields^ 



B t ,n X U,k) 



npxe 



pxfa 



kvr 
Pxh 



when the width bx(j,k) coincides n times the periodicity 
px, as depicted in Fig. [5] (left) for n = 1, taking the 
open orbit A(0, 1) generated by the periodicity pa as an 
example. The positions given by Eq. (|9|) with n = 1 are 
indicated by downward triangles in Fig. [HIb), showing 
that the major features are explicable with the transverse 



FIG. 7. (Color online) Open orbits in the hexagonal lateral 
superlattice. (a) Open orbits in the reciprocal space. Open 
orbits generated by the periodicities p\ , A(0, 1), A(l, 2), and 
by pb, B(0, 1), are shown by thick lines. Black dots and circles 
represent the reciprocal lattice points and the Fermi circles, 
respectively, (b) Illustration in the reciprocal space of the 
index (j, k) of the open orbits, (c) Open orbits in the real 
space. Width b and periodicity along the length c are shown. 



resonances with lower indices (j, k) for the periodicity 
either pa or pb ■ 

Since the open orbits possess the periodicity cx = 
h/(eBpx) (regardless of the indices (j, k)) in the direction 
of the propagation as illustrated in Fig. [TJc) , and the pe- 
riodicity of the potential modulation, p' x , is present also 
in that direction (p' A = a and p' B — y3a), an alternative 
types of the resonances, longitudinal resonances, can be 
considered when cx equals n ■ p' x , namely at 



B, 



npxp'x e ' 



(10) 



as exemplified in Fig. [5] (right). Note that Eq. (fTOj) re- 
duces to 



B, 



(11) 



for both the periodicity pa and pb hi HLSL. A similar 
resonance has been reported in a ID lateral superlat- 
tice slightly modified to have the modulation also along 
the direction perpendicular to the principal axis (strictly 
speaking, therefore, it is a 2D rectangular lateral super- 
lattice), showing a resistivity maximum at the condition 
corresponding to Eq. (IT0|) .— The positions of the reso- 
nances described by Eq. (jlll) are marked by vertical ticks 
in Fig. [Sib), suggesting that some of the minor struc- 
tures in the magnetoresistance are actually originating 
from the longitudinal resonances. Remnant unidentified 
minor structures could possibly be resulting from higher 
order terms (n > 2) in Eq. @ for periodicities pa and 
Pb, or from the transverse resonances for still smaller 



FIG. 8. (Color online) Resonance conditions for the open or- 
bit A(0, 1) generated by the periodicity pa- Resonance takes 
place when the width 6a of the open orbits coincides the pe- 
riodicity pa responsible for the generation of the orbit (trans- 
verse resonance, B — Bt) or when the period ca along the 
length of the orbit matches the periodicity of the modulation 
Pa in the direction of the propagation (longitudinal resonance, 
B = B t ). 



lattice spacings embedded in the hexagonal lattice ex- 
emplified by pc in Fig. EJa). (Positions for longitudinal 
resonances are still described by Eq. (ITTj) for such pe- 
riodicities.) However, these resonances are so densely 
distributed in the magnetic field range shown in Fig. [6] 
that it is rather difficult to unambiguously identify the 
small features in Fig. Hth) with these resonances within 
the resolution of the present experiment. 

Finally, we plot in Fig. |H] the oscillatory part of low- 
field magnetoresistance for a HLSL with a — 100 nm 
analogous to that shown in Fig. Eth) for a = 200 nm. 
We found that both CO and AM of SdHO are extremely 
small for a — 100 nm, hindering us from drawing out re- 
liable information on the potential modulation through 
the analyses similar to those done for a = 200 nm in the 
preceding subsections IIII A l and IIII Bl This is attributed 
to the weakness of the potential modulation, owing to the 
smallness of the period a. Compared with these oscilla- 
tions, GROO can be identified clearer as demonstrated 
in Fig. IH1 This is consistent with the higher sensitiv- 
ity of GROO, compared with CO, to a small periodicity 
p, mainly resulting from a higher characteristic temper- 
ature that governs the decrease of the oscillation ampli- 
tude with increasing temperature or decreasing magnetic 
fields In Fig. [51 two traces taken on different cooling 
downs are shown. Major peaks clearly recognizable in 
both traces take place at the positions, indicated by the 
downward triangles, of the transverse resonances Eq. §§§ 
with n = 1 having lower indices (J, k) for periodicity pa 
or pb, similar to the case in Fig.^b). Slight shifts in the 
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FIG. 9. (Color online) Oscillatory part of the low-field magne- 
toresistance for the HLSL with a = 100 nm obtained by sub- 
tracting slowly-varying polynomial background. Downward 
triangles indicate the positions for transverse resonances given 
by Eq. ((9}. The lower trace is taken with dB/dt — 1 mT/s. 
The upper trace, offset by 0.02 for clarity, is taken on a differ- 
ent cooling down roughly one month after the bottom trace 
was acquired, with 5 times slower sweep rate dB/dt — 0.2 
mT/s thus exhibiting higher resolution. 



peak positions between the two traces are possibly due 
to the small difference in the electron density between 
different cooling downs. Some of the remnant features 
in Fig. [5] are probably resulting from higher n transverse 
resonances or from longitudinal resonances Eq. al- 
though unambiguous identification turned out to be dif- 
ficult as was the case for a = 200 nm. 



IV. DISCUSSION 
A. Towards artificial massless Dirac fermions 

Generation of artificial massless Dirac fermions (MDF) 
in conventional 2DEGs is an enticing possibility. Not 
only does it provide an alternative arena to pursue the 
physics of Dirac fermions, but it also offers an opportu- 
nity to precisely control or modify the properties of Dirac 
fermions; modern semiconductor nano-lithography tech- 
nology, along with the length scale orders of magnitude 
larger than the inter-atomic distance, allows us to manip- 
ulate the hexagonal lattice at will. For instance, we will 
be able to introduce designed strain into the lattice to 
generate an effective magnetic field^ or fabricate nano- 
ribbons truncated at the edges along desired orientations 
(zigzag or armchair), 63 which will be extremely difficult 
to be performed on natural graphcnc. 

Several attempts have been made to generate MDF in 
conventional GaAs/AlGaAs 2DEGs^— reporting (i) nar- 
rowing of a photoluminescence peak^ (ii) small ampli- 
tude magnetoresistance oscillations periodic in B in the 
quantum Hall regime 3 (qualitatively resembling those ob- 



served in the antidot lattices 2 ^), and (iii) splitting of the 
cyclotron resonance absorption. 4 The splitting in (iii) is 
similar to AM of SdHO in the present study in the sense 
that both directly reflect the broadening of the Landau 
levels due to the modulation (Landau bands). Although 
these phenomena undoubtedly stem from the potential 
modulation introduced into 2DEGs, they are not partic- 
ularly sensitive to the lattice type of the 2D modulation. 
Above all, they are not probing the key ingredient in 
the generation of the artificial MDF, the formation of 
minibands by Bragg reflections. We believe that our ob- 
servation of the open orbits engendered by Bragg reflec- 
tions represents one step forward toward the realization 
of MDF. 

However, a number of improvements are still needed 
to be made. (Criteria for achieving MDF in 2DEGs are 
discussed in detail in Ref. First, the electron density 
n e has to be reduced to ~ no = (4/\/3)a~ 2 in order 
to tune the Fermi energy Ep close to the Dirac point. 
The density no equals 6 x 10 13 and 2 x 10 14 m~ 2 for 
a = 200 and 100 nm, respectively. For this purpose, we 
prepared HLSL devices equipped with a backgate which 
allows us to vary n e . 64 However, the values of ap- 
pear to be too small to be achieved with the 2DEG wafer 
used in the present study, preserving the quality (mobil- 
ity) good enough for the minibands to be formed without 
being hampered by the disorder. It will be necessary to 
start with a 2DEG having a higher mobility and a much 
smaller n e . Reducing the lattice constant a will also be 
of help; by reducing a to 50 nm, which is not impossible 
using our high-resolution resist, tid is augmented to a 
more amenable value of 9 x 10 14 m~ 2 , 6 - Second, in order 
to avoid the Dirac cones from energetically overlapping 
with other minibands, the sign of the hexagonal potential 
modulation ought to be positive (repulsive) Unfor- 
tunately, none of the oscillations discussed in the present 
study (CO, AM of SdHO, and GROO) is sensitive to the 
sign of the modulation, and we therefore do not know 
whether our hexagonal potential modulation is repulsive 
or attractive. If the strain induced by the resist turns out 
to introduce attractive potential, it will be necessary to 
switch to the honeycomb lattice modulation dual to the 
hexagonal latticed 

Very recently, Gomes et al. reported^ the generation 
of MDF on a 2DEG at a copper surface, using carbon 
monoxide molecules as a negative gate working on the 
2DEG. The molecules were assembled by the atomic ma- 
nipulation technique employing a scanning tunneling mi- 
croscope (STM). They probed the resulting Dirac cones 
using the scanning tunneling spectroscopy. Although 
their result is impressive, their method requires acrobatic 
operation of the STM. Also, it will be probably not very 
easy to perform transport measurement on a 2DEG at 
the copper surface residing just above the bulk of the 
copper. We therefore believe that generating MDF in 
a conventional semiconductor 2DEG is a challenge still 
worth pursuing. 
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B. Towards the observation of the bubble phase 

We extended the magnetoresistance measurement in 
the dilution refrigerator (15 mK) up to 9 T, in search of 
the effect induced by the hexagonal modulation in the 
quantum Hall regime. We especially focused on 1/4 and 
3/4 fillings of the N = 2 and higher Landau levels, seek- 
ing for signals related to the bubble phase theoretically 
predicted to be the ground state at these fillings^— 

The bubble phase is a charge density wave (CDW) 
state in which clusters of two or three electrons are ar- 
ranged in the hexagonal lattice. Experimentally, reen- 
trant integer quantum Hall effect (RIQHE) observed at 
filling factors v = 17/4, 19/4,... was interpreted as re- 
sulting from the bubble state pinned by impurities.^ 
Later, microwave resonances were observed at these fill- 
ings, which were ascribed to the pinning mode of the 
bubble stated For these experiments, 2DEGs having a 
very high mobility (/x > 1000 m V s ) were required, 
since the formation of the fragile bubble state is readily 
quenched by disorder. On the other hand, external mod- 
ulation having the same lattice constant and the lattice 
type, namely the hexagonal lattice, as the bubble state 
is theoretically predicted to stabilize the bubble phased 
We therefore expect the bubble states to be formed in 
optimally designed HLSLs even if the mobility is not as 
high. In fact, anisotropic magnetotransport was observed 
to be induced by external modulation in ID lateral su- 
perlattices with a mobility ^ ~ 100 m 2 V _1 s _1 , which 
was ascribed to the stabilization of the stripe phase ; 1 ^ 11 
similarly fragile phase predicted to be the ground state 
at the half fillings. 7-9 Stabilization of the bubble phase 
by the external modulation in HLSLs would provide di- 
rect information on the spatial distribution of the charge 
through the lateral size and the shape of the introduced 
modulation. Note that experimental findings related to 
the bubble phase reported so fa r 66 ' 67 do not have sensi- 
tivity to the spatial distribution. 

Unfortunately, we have found no changes above 1 T 
attributable to the introduction of the modulation in 
HLSLs thus far for both a = 200 and 100 nm. This is 
probably because either the introduced lattice constant 
a or the strength of the potential modulation was not 
appropriate. The lattice constant ocdw of the bubble 
state is theoretically predicted to be acDW ~ 3i? c = 
ll. 7 ' 8 For our n e = 3.9 x 10 15 m" 2 , a CD w = 88, 
93, 98, and 102 nm for v = 17/4, 19/4, 21/4, and 23/4, 
respectively. Therefore, a = 200 nm appears to be rather 
too large to promote the formation of the bubble phase. 



On the other hand, a = 100 nm roughly matches the 
predicted ocdw- I n this case, however, the modulation 
was probably too weak to overcome the detrimental effect 
exerted by the disorder. 

Since the magnetic length I hence the lattice constant 
Ocdw at a fixed filling factor v = n e h/(eB) increases 
with decreasing n e , we can increase the suitable a hence 
the modulation amplitude by using smaller n e . There- 
fore, it will be desirable, here again, to prepare HLSL 
samples with 2DEGs with higher mobility and smaller 



V. CONCLUSIONS 

We have observed three types of oscillatory phenom- 
ena in the magnetoresistance of hexagonal lateral super- 
lattices (HLSLs): commensurability oscillations (CO), 
amplitude modulation (AM) of Shubnikov-de Haas os- 
cillations (SdHO), and the geometric resonances of open 
orbits (GROO). Both CO and GROO contain compo- 
nents deriving from two periodicities, pa = Via/2 and 
Pb = <i/2, immanent in the hexagonal lattice with the lat- 
tice constant a, while only the larger periodicity pa man- 
ifests itself in the AM of SdHO. As in the case of square 
or rectangular 2D lateral superlattices, amplitude of CO 
in HLSL is much smaller than that in ID lateral superlat- 
tices having similar strength of the potential modulation. 
By contrast, magnitude of AM in the SdHO is compara- 
ble to that in ID lateral superlattices. The observation of 
GROO reveals that minibands are generated by Bragg re- 
flections from the hexagonal superlattices, which requires 
highly ordered modulation potential with a small enough 
lattice constant close to the Fermi wavelength. With fur- 
ther adjustment of the parameters of the 2DEG and of 
the modulation potential, it might become possible in the 
future to probe by magnetotransport experiments artifi- 
cially designed massless Dirac fermions (MDF) generated 
in the miniband structure. 
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